Possible effects of charge frustration in Na^CoC^: bandwidth suppression, charge 
orders and resurrected RVB superconductivity 
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Charge frustration due to further neighbor Coulomb repulsion can have dramatic effects on the 
electronic properties of Na^CoC^ in the full doping range. It can significantly reduce the effective 
mobility of the charge carriers, leading to a low degeneracy temperature tp < T. Such strongly 
renormalized Fermi liquid has rather unusual properties — from the point of view of the ordinary 
metals with ef > T — but similar to the properties that are actually observed in the Na^CoC^ sys- 
tem. For example, we show that the anomalous thermopower and Hall effect observed in Nao.7Co02 
may be interpreted along these lines. If the repulsion is strong, it can also lead to charge order; 
nevertheless, away from the commensurate dopings, the configurational constraints allow some mo- 
bility for the charge carriers, i.e., there remains some "metallic" component. Finally, the particularly 
strong bandwidth suppression around the commensurate x = 1/3 can help resurrect the RVB su- 
perconductivity, which would otherwise not be expected near this high doping. These suggestions 
are demonstrated specifically for a tJ-like model with an additional nearest neighbor repulsion. 



I. INTRODUCTION 

Recent discovery^ and confirmation 2 ' 3 ' 4 of supercon- 
ductivity in Nao.35Co02 • I.3H2O has stimulated many 
studies of this material, and in particular of its un- 
hydrated precursor host material Nao.7Co02- The 
Na^CoC^ series has been known for more than five 
years 5 ' 6 ' ' - 8 for its unusual transport properties such as 
large thermoelectric power and nearly linear- T depen- 
dence of the resistivity, indicative of strong correlation 
effects. This has been brought up even further by a num- 
ber of recent careful experiments piLifiiiAiiSiiiiii 

It appears from these studies^ that there is a large 
suppression of the valence band width (equivalently, large 
effective mass enhancement) — by an order of magnitude 
compared with the LDA band structure calculations^ 

We suggest that such large renormalization may be 
caused by strong Coulomb repulsion between charge car- 
riers on neighboring Co sites, and point out that a num- 
ber of unusual properties of the system may be explained 
by the resulting low fermion degeneracy temperature. We 
also consider other possible effects of such repulsion, in 
particular, charge ordering. 

The plan of the paper is as follows. To be specific, we 
consider a tJ-like model with additional strong nearest 
neighbor repulsion V. In sections IlllVI we concentrate 
on the dominant t, V energetics. The study is done by 
considering Gutzwillcr-like trial fcrmionic wavefunctions 
(projected Fermi liquid) with additional nearest neighbor 
correlations input through a Jastrow-type configurational 
weighting factor. The strength of the input correlations 
serves as a variational parameter. 

Sec. IIIII studies the properties of these wavefunctions. 
We find that up to moderate input correlations, the wave- 
function indeed describes a renormalized Fermi liquid, 
consistent with the initial motivation. We also realize 
that for strong input correlation and over the doping 
range 0.27 < x < 0.5 and 0.5 < x < 0.73, our Jastrow- 



Gutzwiller wavefunction has a a/3 x a/3 charge order, 
which is inherited from the charge distribution proper- 
ties of the classical Jastrow weight on the lattice. Since 
such a state is beyond our initial motivation, we examine 
its properties and treat it very critically whenever the 
variational parameter is driven into this regime. 

In Sec. lIVI we develop a convenient renormalized mean- 
field picture for the energetics in the entire doping range. 
We identify the regime of the renormalized Fermi liq- 
uid state, and also the regime where the strong repulsion 
drives the optimal Jastrow-Gutzwiller wavefunction into 
the -\/3 x a/3 charge order state. 

In Sec. El we confirm the renormalized meanfield pic- 
ture with numerically accurate evaluations with the trial 
wavefunctions. We also perform a more detailed study 
of the possible a/3 x a/3 charge order by comparing with 
the more conventional CDW states. We find that our 
Jastrow-Gutzwiller wavefunction in the a/3 x a/3 regime 
is rather good energetically and suggest some ways for 
improving the energetics further. 

In Sec. IVII we add the J term and consider the is- 
sue of RVB superconductivity at low dopings. This is 
done with the help of the renormalized meanfield pic- 
ture. Without the Jastrow renormalizations, the RVB 
superconductivity would not survive to the experimen- 
tally observed x = 0.35. We find that the bandwidth 
suppression due to charge frustration may indeed res- 
urrect the superconductivity near x = 1/3 where such 
renormalizations are strongest, particularly if we allow 
the coexisting charge order. We speculate that this may 
be relevant to explain the narrow doping range in which 
the superconductivity has been found^ 

Finally, in Sec. lVIll we conclude with some simple pre- 
dictions for the experiments from the developed charge 
frustration picture. Most notably, transport properties 
such as thermopower and Hall effect of the Fermi liquid 
with low degeneracy temperature resemble those of the 
Nao.7Co02 system; these properties look rather unusual 
from the perspective of conventional metals. 
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Before proceeding, we remark about the possibility of 
charge order&i 6 . in Na^CoC^ . The experimental situa- 
tion is not settled on this issuei 12 i 13 ' 14 We favor the pic- 
ture where there is no charge order, but only strong local 
correlations. Charge ordering transition should exhibit 
itself in an abrupt change in transport properties, which 
has not been observed. 

In the present paper, we do spend a lot of time dis- 
cussing the particular vo x \/3 order, since it inevitably 
arises in our systematic treatment of the concrete model. 
We should warn the reader that the details are likely 
strongly model-dependent. Since we do not know the 
precise microscopic model, the presented analysis of the 
charge order should be viewed only as an initial sketch of 
what might happen. The reported work is done with the 
nearest neighbor repulsion (and nearest neighbor Jastrow 
correlation) only. Including further neighbor correlations 
would frustrate the x charge order and extend the 
renormalized Fermi liquid regime, but might also lead to 
more complicated charge orders. We have not pursued 
such studies systematically, concentrating on the nearest 
neighbor case only. 



II. tV MODEL. JASTROW-GUTZWILLER * 

For concreteness, we consider the following single 
band Hubbard model with additional nearest neighbor 
repulsion 16 on a triangular lattice 



H, 



LV 



(ij) 
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+ h.c.)P G + V^n i n j 

(ij) 



(1) 



Large onsite repulsion is taken into account using the 
Gutzwiller projector P G to project out double-occupation 
of sites — this gives the "f part" as in the familiar t J 
Hamiltonian. We focus primarily on the effect of adding 
strong nearest neighbor charge repulsion V, and refer 
to Eq. as tV Hamiltonian. (We will consider the 
full tJV Hamiltonian with J <C t,V later.) The band 
is less than half-filled, with the average fermion density 

(4t C «t) = P < 1. 

In the context of the Na^CoC^ system, c| CT creates a 
spinful hole and represents the motion of a Co 4+ (S = 
1/2) site, while Co 3+ (5 = 0) sites have no holes; p — 
1 — x. This is shown schematically in Fig. ^ F° r more 
details, see also Refs. Ilmlall9l We take t > 0, as in 
Refs. TT^fTgL This is consistent with the photoemission 
studiesfii and also with the high temperature behavior 
of the Hall coefficien1p u i (see Sec. lVII|) . In this picture the 
end compound NaCo02 consists of all Co 3+ with no holes 
(p = Q,x = 1), whereas the hypothetical end compound 
Naa; = oCo02 consists of a Mott insulator with all Co 4+ 
sites each carrying S = | (p = 1, x = Q). From the latter 
point of view, Na^OoC^ can be viewed as electron doping 
by a concentration of x electrons into a Mott insulator. 

Here, in order to gauge the effect of the nearest neigh- 
bor repulsion, we perform a trial wavefunction study of 



Co 4+ Co 



,4+ 



3+ 




FIG. 1: Na I Co02 : Schematic pictures (borrowed from 
Ref. explaining the single band electronic model. The 
charge carriers are spin- 1/2 charge q — ej holes of density 
p = 1 — x. 



this strongly correlated system. When V — 0, a good 
trial wavefunction is obtained by Gutzwiller-projecting a 
simple free fermion state: 

|*g) = Pg\*o) = det c fliT ■ ■ • Cfl N/2 T 

{R},{R'} 



det [<pi{R'j)] cg^i ■■■c R 



N/2- 1 - ' 

where the sum is over all configurations of spin up and 
down fermions with all Rj and R'j, distinct, {R}C\{R'} = 
0. Away from half-filling, this Gutzwiller wavefunction 
is a Fermi liquid state; this can be confirmed, e.g., by 
measuring the quasiparticle Z from the step in (fife). In 
this state, we have approximately {riirij) ~ p 2 , while the 
fermion kinetic energy can be fairly accurately estimated 
from that of the preprojected free fermions22i2L2S 



(y G \H t \^ G ) 
(*g|*g) 



9f 



o 



with 



9t 



1 



2.T 



l-p/2 1 



(2) 



(3) 



This is commonly referred to as the Gutzwiller approx- 
imation, or the renormalized mean field theory. The 
renormalization factor g t can be obtained by counting the 
number of real space configurations available for hopping 
in the projected and preprojected states, and ignoring all 
other wavefunction differences. 

Turning on the nearest neighbor repulsion V, we 
schematize its effect on the ground state by introducing 
an additional Jastrow-type factor 



exp 



W 



(4) 



for each real space configuration of fermions in the above 
Gutzwiller wavefunction. W > effectively suppresses 
the nearest neighbor occupation probability, and can be 
varied to optimize the trial energy of the tV Hamiltonian. 
We will refer to this wavefunction as Jastrow-Gutzwillcr 
(JG) <F /G . 
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It is clear that the effect of V is most severe for x = 1/3 
and 2/3. For V ^> t we expect W ^S> 1, in which case the 
vacancies (for x = 1 /3) or the spin carrying holes (for x = 
2/3) would form a structure to minimize the 

repulsion. Furthermore in this state the particles cannot 
hop without paying the energy V. For intermediate V 
and away from commensuration some remnant of this 
"jamming" phenomenon may remain and this is what we 
would like to investigate in this paper. 



III. PROPERTIES OF * /G . CLASSICAL 
LATTICE GAS SYSTEM. 

Before presenting the optimized energetics with tyjc, 
we first discuss its properties. For nonzero but small W, 
it still describes a Fermi liquid state with some further 
rcnormalizations compared to the Gutzwillcr wavefunc- 
tion. However, one has to be careful when W becomes 
large: The probability of finding a particular configura- 
tion of charges now has an additional "classical" weight 



exp[-W/£ 



(ij) 'H"jJ, 



and one has to be wary of the possi- 



bility of phase transitions in the corresponding statistical 
system. Factoring out many possible spin assignments for 
each charge configuration, we need to consider a classical 
system of particles with nearest neighbor repulsion on a 
triangular lattice with the classical partition function 



^class = ]T e"^ 



{«»} 



(5) 



Here rii = 0, 1, and we work at fixed density p as appro- 
priate for the discussion of our trial wavefunctions with 
fixed fermion number. W plays the role of the inverse 
temperature in this classical system, W = T _1 (the clas- 
sical repulsion strength is set to one). 

This lattice gas system has been extensively studied in 
statistical physicsj2&2i most notably as a model for ad- 
sorbed monolayers of rare-gas atoms on graphite. Also, it 
is equivalent to a triangular lattice Ising antiferromagnet 
in an external field; fixed particle density corresponds to 
fixed magnetization in the Ising system. 

The phase diagram in the p — T plane is shown in 
Fig. |21 It is symmetric with respect to p — 0.5 due to 
particle-hole symmetry in this classical system, and we 
discuss the p < 0.5 part only. 

At high temperatures (small W) the system is in a 
disordered gaseous phase. For small particle density 
p < 0.27 the system remains in the gas phase all the 
way to zero temperature. For densities 0.27 < p < 0.5 
the system "crystallizes" into a V3 x \/3 state at low 
temperature. This state is characterized by a preferen- 
tial particle occupation of one of the three sublattices 
of the triangular lattice; the order is strongest near the 
commensurate p = 1/3, 2/3. Note that away from these 
p the \/?> x \/?> phase is more properly characterized as a 
density wave state rather than a crystal. In particular, a 
fraction of particles remains relatively mobile having no 
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FIG. 2: Phase diagram of the classical lattice gas system 
Eq. (Jy . Note the symmetry relative to p — 0.5 due to particle- 
hole symmetry. 



activation energy for their motion. The W = oo model 
is equivalent to the Baxter's hard hexagon model and is 
exactly solvable i^i 

It is useful to have the following caricatures of the 
charge motion in the \/3 x y/3 phase, appropriate at 
low temperature (large W) and near the commensurate 
filling 1/3. For p > 1/3, we have one sublattice, say 
A, completely occupied by particles (these particles are 
almost localized), while the remaining small density is 
smeared out relatively uniformly over the honeycomb lat- 
tice formed by the B and C sublattices. For p < 1/3, we 
picture the B and C sites as completely empty (with 
almost no density fluctuation), while all particles are 
spread over the A sublattice. The charge motion is 
achieved by hops from occupied A sites to neighboring 
empty A sites via B or C sites; the most effective such 
hops involve at least two neighboring empty A sites in 
order to avoid the repulsion energy cost for the interme- 
diate step (see Fig.|S]in Sec. Ivj) . 

Returning to our trial wavefunction, we expect it to 
roughly inherit the charge distribution properties of the 
lattice gas system Eq. |JS}. Thus, for W such that the 
classical system is in its disordered phase, the wavefunc- 
tion realizes a Fermi liquid state. On the other hand, 
when the classical system is in the x phase, we 
have checked by variational Monte Carlo studies that the 
resulting wavefunction has a charge density wave (CDW) 
order but retains some liquid properties. (The transi- 
tion point W c is slightly different for the wavefunction 
and the classical system.) Note that the above Jastrow 
weight realizes a soft projection satisfying the nearest 
neighbor repulsion V. It follows that a hard such pro- 
jection (W — ► oo) leads to the y3 x \/3 density wave for 
0.27 < p < 0.5 and 0.5 < p < 0.73. 

Before proceeding further, we want to emphasize again 
that we set out to study Fermi liquid renormalizations 
induced by nearest neighbor repulsion. Fermi liquid state 
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is achieved in the parameter regime where the classical 
system remains disordered; in this case, our treatment is 
consistent. 

On the other hand, when the nearest neighbor repul- 
sion is strong, it drives the optimal variational param- 
eter W of \P jg into the regime with the CDW order. 
In this case, we have to be very cautious in interpreting 
the "transition" and the resulting state, since our origi- 
nal assumptions about the properties of the wavefunction 
no longer hold. We may still interpret this as a sign of 
an instability towards a different state (most likely with 
charge order), but the JG wavefunction in this regime 
should be treated very critically, particularly since it has 
somewhat unusual charge distribution properties. Thus, 
one should at least examine other more conventional trial 
states with different orders. This is done in Sec. |Vj 



IV. ENERGETICS WITH * JG : 
RENORMALIZED MEANFIELD PICTURE. 

We now proceed to the actual energetics with the 
Jastrow-Gutzwiller trial wavefunctions. It is possible 
to perform essentially exact evaluations of the expec- 
tation values with such fermionic wavefunctions using 
a well established and documented Variational Monte 



Carlo (VMC) procedure. 21 ! 22 ! 25 Such detailed studies of 
the tV energetics at experimentally relevant x — 0.7 and 
0.35 are reported in the next section. 

It is also possible to study more complicated Hamilto- 
nians. However, VMC evaluations are computationally 
rather costly. Furthermore, they become inconclusive 
when the energy differences become very small. This 
is particularly the case when we attempt to study the 
physics at energy scales below the dominant t and V, 
e.g., if we want to resolve the spin sector, or study pair- 
ing instabilities due to the J term. 

Useful and fairly accurate guidance is obtained through 
the following "renormalized meanfield" procedurerL 2 ^ 2 ^ 
which is much simpler computationally and is also more 
amenable to interpretation and extrapolation in the 
regime where VMC results become inconclusive. Gen- 
eralization of the configuration counting arguments men- 
tioned earlier leads to the following estimate of the hop- 
ping energy renormalization in the JG wavefunction rel- 
ative to the unprojected free fermion wavefunction: 



Wj G |cT T c 3 - T |Wj G ) 

(*,7G|*JG) 



^o|cIt c jtI^o) 
(*o|*o) 



(6) 



with 



J 



9t[hJ\ 



p(l - p/2) 



S(rii — 0)S(rij — 1) exp 



^ [Ucl ass [n i = l,n j =0] -Ucl ass [ni=0,7lj = l] j 



r 



(7) 



Here, ((. . .)) denotes averaging in the classical lattice 
gas system with the weight ~ exp(— ll c i ass [n]) discussed 
earlier. When obtaining this expression, similar to the 
original Gutzwiller approximation Eq. J3J), we again ig- 
nored the details of the fermionic determinant weight- 
ing of configurations, but kept the Jastrow weighting. 
Only configurations with the occupied j site and un- 
occupied i site contribute, and the specific "transition 
weight" comes from the corresponding Jastrow weight- 
ing of the configurations before and after the hop. Note 
that il c i ass [nj = l,rij = 0] - it c i ass [n l = 0, n 3 = 1] is a lo- 
cal energy term involving only the affected sites i, j, and 
their immediate neighbors. 

Similarly, we can approximate the nearest neighbor re- 
pulsion energy by 



jG\riinj\y jg 



JG 



(8) 



The required classical expectation values are readily eval- 
uated via a Monte Carlo study of the lattice gas sys- 
tem. As we will see in the next section, such renormal- 
ized meanfield procedure indeed gives fairly accurate esti- 
mates of the expectation values in the Jastrow-Gutzwiller 



wavefunction, both in the metallic and the density wave 
regimes. 

We can now develop an overall picture for all fermion 
densities. Particular cuts through the results are shown 
in Figs. 13 and 0] Figure [3] shows the hopping renormal- 
ization factor g t as a function of x for a number of fixed 
W. (Complimentary cuts through the data for the spe- 
cific fixed x = 0.70 and x = 0.35 can be also found in the 
next section.) Figure is the core of the present paper. 

The W = curve gives precisely the original Gutzwiller 
approximation Eq. for the no-double-occupancy con- 
straint. This sets a useful reference for gauging the addi- 
tional effect of the nearest neighbor repulsion. The curve 
with the largest W = 8, on the other hand, essentially 
realizes a complete projection that satisfies the nearest 
neighbor repulsion; this is the maximal renormalization 
that can be achieved with such nearest neighbor corre- 
lations. The phase boundary of the classical lattice gas 
(cf. Fig. [2J is sketched by a thick dark line: All points 
above the line are in the disordered phase (Fermi liquid 
wavefunctions), while points below the line are in the 
V3 x \/3 density wave phase. 

Figure 01 shows a similar plot for the repulsion en- 
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Renormalization factor for hopping 




tV 'phase diagram' with *¥jq in renorm. meanfield 



0.4 0.5 0.6 
doping X 

FIG. 3: Jastrow-Gutzwiller renormalization factor for hop- 
ping, Eq. 0, as a function of doping for a number of fixed 
W. Evaluations are done via classical Monte Carlo study of 
the lattice gas system. The dark thick line delineates the 
phases of ^ jg- The y/S x y/3 phase lies below the thick line 
and in this region the exhibited gt is averaged over all bonds. 




0.2 0.3 0.4 0.5 0.6 
doping X 



1.0 



FIG. 5: "Phase diagram" obtained by optimizing the tV 
Hamiltonian over the JG wavefunctions. The different 
"phases" correspond to the difference in the physical prop- 
erties of <&jg as a function of W . Calculations are performed 
in the renormalized meanfield approximation using the data of 
Figs.ElBl an d t ne f ree fermion (H t )o(x). The phase diagram 
is not expected to be symmetric relative to V — cf. Fig. [3] 
The observed rough symmetry is due to compensating ten- 
dencies in gt(x) and (Ht)o(x) that make the actual kinetic 
energy "more symmetric" relative to x = 0.5. 



E nn [W] - E nn [W^ co] 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



doping X 

FIG. 4: Nearest-neighbor repulsion energy per site, referenced 
to the W — > oo value. The y/3 x \/3 charge ordered phases of 
^ jg lie below the thick line. 

ergy E nn = Y^(ij)(( n i n j)) P er site ( cf - Eq.|SJ), which we 
reference to the minimal possible repulsion energy at a 
given density: E nTlt w^oo — for p S [0, 1/3], ip — 1 for 
p € [1/3, 2/3], and 6p - 3 for p € [2/3, lj. When plot- 
ted in this way, the result is symmetric with respect to 
p = 0.5 due to classical particle-hole symmetry. Again, 
the classical phase boundary is sketched with a thick dark 
line. Observe that the curves with W > 5 give almost 
complete "minimum-nearest-neighbor" projection. 

With this data, and also using the free-fermion 
(Ht)a(x) (not shown), we can optimize the full tV Hamil- 
tonian in this renormalized meanfield procedure for ^ jq- 
The resulting "phase diagram" can be seen in Fig. [S] 



For each doping x 6 [0.27, 0.73] we show the "critical" 
V/t that drives the optimal W into the regime with the 
V3 x V3 order. 

We emphasize that this "phase diagram" is for the 
optimized Jastrow-Gutzwiller wavefunction only; in par- 
ticular, the exhibited "phase transition" corresponds to 
the transition in the properties of ^ jg as a function 
of parameter W. It need not correspond to the actual 
phase diagram of the tV Hamiltonian. It is exhibited 
here primarily to delineate the regimes where the JG 
renormalized Fermi liquid can adequately describe the 
tV model energetics, and also where such description is 
no longer possible. In the latter case, one should se- 
riously examine other physical states paying particular 
attention to charge order. Whether the JG wavefunction 
in its y/3 x phase can adequately describe the possi- 
ble charge ordering in the system is a separate question 
that requires a detailed study. We will discuss this more 
specifically in the next section. Here we only note that 
it is rather fortuitous that our trial wavefunction with 
a single variational parameter exhibits two phases, and 
the initially "unexpected" charge-ordered state should be 
treated with great caution. 

We now return to the main question of this work — the 
bandwidth suppression due to nearest neighbor repulsion. 
Again, consider Fig. El A conservative approach is to in- 
sist that we consider Fermi liquid wavefunctions only. In 
this case, we should disregard the data points that end up 
in the charge ordered phase. We still see that there can 
be significant renormalizations by a factor of three to five 
relative to the bare hopping amplitude even remaining in 
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the Fermi liquid state. For a fixed W, these renormaliza- 
tions are strongest near the commensurate 1/3 and 2/3 
fillings, and weakest near x = 1/2. Also, as can be im- 
plied from the "phase diagram" , the effect of the nearest 
neighbor repulsion V is strongest near x — 1/3, 2/3. 

On the other hand, if we are to take the Jastrow- 
Gutzwiller y/3 x \/3 density wave regime seriously, there 
can be even stronger renormalizations of the hopping en- 
ergy, particularly near the commensurate densities. As 
we will suggest in the following more specific discus- 
sion of the CDW regime, the entire picture provided by 
Fig. 13 including the data under the phase boundary is 
indeed useful, but may require some less important ad- 
justments. This is because in the V3 x regime the 
charge order is such that there remain mobile (even if 
strongly constrained) carriers; there is no charge gap 
since there is no nesting for the considered dopings. The 
Jastrow-Gutzwiller wavefunction and the above renor- 
malized meanfield treatment also capture this, and give 
a first useful guess on the effect of charge order on the 
fermion kinetic energy. 
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FIG. 6: Expectation values of the hopping and nearest neigh- 
bor repulsion energies in the Jastrow-Gutzwiller wavefunc- 
tion at doping x — 0.70 (fermion density p — 0.3). We also 
show the corresponding renormalized meanfield values. Verti- 
cal arrow near W c « 3.3 indicates the transition point in the 
classical lattice gas system. 



V. ENERGETICS WITH * JG : VMC STUDY. 
POSSIBLE CDW. 

We now consider the tV model energetics in more de- 
tail for specific x = 0.70 and x — 0.35. These values 
are relevant for the unhydrated and hydrated Na x Co02. 
The evaluations with the trial wavefunctions are done 
essentially exactly using VMCiSSiSit This more concrete 
setting will allow us to discuss some robust features that 
emerge from our study vs the specifics of the particular 
Hamiltonian used to model charge frustration. Since an 
accurate treatment of the CDW states may depend on 
specific details, the present discussion is only intended to 
give a flavor of the possibilities that should be considered. 



A. Doping x = 0.70 

Figure shows expectation values of the two parts of 
the tV Hamiltonian in the Jastrow-Gutzwiller wavefunc- 
tion for varying W evaluated using VMC. It also shows 
the renormalized meanfield approximation to these ex- 
pectation values. As mentioned in the previous section, 
this approximation is indeed fairly accurate and can be 
taken seriously. Since we will be comparing several trial 
states, we will use only VMC results in this section. 

At this particle density we have W c « 3.3 in the cor- 
responding lattice gas system. Notice that near W c the 
repulsion energy drops quickly and essentially all the way 
to zero, similar to the transition in the classical system. 
This is because it is possible to completely satisfy the 
repulsion energy by arranging charges so that there are 
no nearest neighbors. Also, such arrangements still allow 
some fermion hopping, so there remains nonzero kinetic 
energy gain even for very large W . 



Optimized tV energetics 




2 4 6 8 10 12 



FIG. 7: Optimization of the tV Hamiltonian over Jastrow- 
Gutzwiller and conventional CDW trial wavefunctions. Inset 
shows the optimal W op t for the Jastrow-Gutzwiller wavefunc- 
tion; arrow indicates the critical W c « 3.3 in the correspond- 
ing classical lattice gas system. 

Using the above data, we can optimize the total en- 
ergy for different values of V/t. The result is indicated 
in Fig. [7| For V/t < 6.0, the optimal W opt remains 
< 2.0 and the wavefunction is metallic with relatively 
weak renormalizations. For larger V/t, the optimal W op t 
jumps to the \/3 x a/3 ordered side, and ^jg has the 
corresponding density wave order. Note that the opti- 
mal W op t remains fairly close to the critical value. In 
this way, while the repulsion energy is almost completely 
satisfied, the system still gains from some of the original 
kinetic energy. 

We now discuss the regime of large V, where the 
Jastrow-Gutzwiller energetics suggests charge ordering. 
First, it is instructive to compare this with the energet- 
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ics in a more conventional CDW trial state. Such a state 
is obtained, for example, by considering a CDW mean- 
field Hamiltonian 

H clw = -^2{ t ij c \<j c j<y +h.c.) -^2A Q cos(Q • n)hi . 

(ij) 1 

(9) 

Aq = Acdw(Q) is a CDW order parameter at the order- 
ing wavevector Q. Here, Acdw serves as a variational 
parameter for the trial wavefunction. In the -\/3 x y3 
phase, the trial -ffcDW has onsite potential — 2A on the 
preferred A sublattice and +A on the B and C sublat- 
tices of the triangular lattice. The trial wavefunction 
is obtained by Gutzwiller projection of the meanfield 
ground state. 

The optimized tV energetics for such more conven- 
tional CDW wavefunction is also shown in Fig. [7| For 
V/t < 4.5 the optimal wavefunction has Acdw ~ 0, but 
develops strong CDW order for larger V/t. In this con- 
ventional CDW state at this filling, we have a coexistence 
of the charge order and Fermi liquid. 

From Fig. [7| we see that the JG wavefunction per- 
forms significantly better than the conventional CDW 
wavefunction. This is simple for the metallic side, since 
the JG wavefunction has an additional variational pa- 
rameter to optimize local correlations compared to the 
plain metallic state with Acdw = 0. On the other hand, 
on the charge-ordered side the Jastrow-Gutzwiller wave- 
function performs better almost entirely due to better 
kinetic energy. As discussed earlier, the JG state retains 
some of the metallic kinetic energy even in the large W 
limit. At the same time, the conventional CDW wave- 
function localizes the fermions to the A sublattice very 
strongly and loses essentially all kinetic energy: in the 
limit of large V , the optimal Acdw ~ V and the opti- 
mal total energy is ~ — t 2 /V. Even though the lowest 
band remains only partially filled, its bandwidth goes to 
zero in the limit of large Acdw • 

Thus, we conclude that the Jastrow-Gutzwiller wave- 
function with the y/3 x V3 order performs fairly well for 
large V. However, this is by no means the end of the story 
even for the tV model. The most serious reservation here 
is that we have not explored other competing states in 
the system for large V. We will not address this. We 
still hope that our approach captures the relevant local 
energetics in the system. 

In the present context, we can explore the energetics of 
the v3 x v3 ordering more systematically. As discussed, 
the complete minimum-nearest-neighbor projection leads 
to the Vo x Vo order. For p < 1/3, we essentially have 
charges living on the A sublattice only and moving pri- 
marily via A — B — A or A — C — A routes, while the 
bonds B — C are rarely used (see Fig. |SJ. In the above, 
we were projecting the uniform free fcrmion triangular 
lattice hopping ground state, while it is clear that in 
the resulting charge ordered state the hops B — C are 
poorly utilized, and more generally the kinetic energy — 
the driving force for uniformity — is less important. In 
the ^/i x \/3 regime, it then seems more appropriate to 
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FIG. 8: a) Schematics of the y/S x \/3 charge order for p < 
1/3. Charges occupy the A sublattice and spend very little 
time on the B and C sublattices. The remaining empty A sites 
can be utilized for charge hopping. There is an intermediate 
repulsion energy cost of V to move an isolated such site, but 
no such cost for two neighboring empty A-sites as shown in 
the figure, b) Dice lattice hopping ansatz motivated by the 
observation that hops B-C are rarely used. 



project a hopping state with strong A — B and A — C 
hopping amplitudes and weak B — C hops. The lim- 
iting case is the dice lattice hopping shown in Fig. |SJd; 
the six-coordinated sites are the A sites, while the three- 
coordinated sites are the B and C sites. The dice hopping 
state by construction has \/3 x -\/3 order. It is easy to 
verify that in the lowest band one half of the fermion den- 
sity is located on the A sublattice. Clearly, this has bet- 
ter repulsion energy than the triangular hopping state. 
To further optimize the nearest neighbor correlations, we 
can introduce a Jastrow weighting as for the uniform hop- 
ping state (note that in the dice case the classical lattice 
gas transition is no longer relevant since the charge sys- 
tem has the \/3 x -\/3 order from the outset). The opti- 
mized tV energetics is shown in Fig.EI and we indeed find 
that the dice hopping ansatz is somewhat better than the 
uniform state. 

Finally, we should point out that we have completely 
ignored the spin physics by considering only unpolar- 
ized wavefunctions. It should be clear that since the 
bandwidth becomes so narrow, there will be significant 
degeneracy — on the tV energy scale — in the spin sector. 
This degeneracy will be resolved in some way or other 
at lower energy scale, and the details will depend largely 
on the specifics of the microscopic Hamiltonian. As an 
example, trying out spin-polarized Jastrow-Gutzwiller 
wavefunctions in the tV Hamiltonian, we find that in 
the charge-ordered regime the fully polarized wavefunc- 
tion performs only slightly worse than the unpolarized 
one. For the dice hopping ansatz, on the other hand, 
the spin-polarized wavefunction performs better than the 
unpolarized one. One can get some feeling of the slight 
differences by examining Fig.|5] Such itinerant ferromag- 
net tendencies become even more pronounced at lower 
fermion density (higher x). 

After the presented detail, it should be clear that the 
energetics can be rather subtle and model dependent, 
particularly in the charge order regime. We now want 
to separate out which features are more robust than the 
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FIG. 9: This is a blow-up of Fig. focusing on the \/3 X V3 
regime and showing additional Jastrow-Gutzwiller type trial 
wavefunctions for the tV Hamiltonian. Besides the unpolar- 
ized triangular lattice hopping ansatz, we also show the opti- 
mized energetics for the corresponding fully polarized state, 
and also for the dice lattice hopping ansatz. 



above specifics. This is important since at present we do 
not have a good knowledge of the microscopic Hamilto- 
nian in the Na^CoC^ system. 

First of all, we conclude that there can be significant 
renormalizations in the metallic wavefunction. The hop- 
ping can be effectively suppressed roughly by a factor 
of three (see Fig. EJ), with the wavefunction retaining 
its Fermi liquid character. The achievable Fermi liquid 
renormalizations may be even larger if we include further 
neighbor repulsion, since this will frustrate the a/3 x a/3 
charge order and give more parameter space to the liq- 
uid state with uniform charge distribution. As long as 
the system remains uniform, this is not sensitive to the 
microscopies. (At this local "high energy" level of anal- 
ysis we completely disregard the low energy instabilities 
of the resulting Fermi liquid state.) 

Our second observation is about the nature of possi- 
ble charge orders in such strongly frustrated system. Our 
JG wavefunctions offer an interesting possibility of essen- 
tially satisfying the nearest neighbor Coulomb repulsion 
while retaining some kinetic energy gain and metallic- 
ity. Projecting the triangular or dice lattice hopping is 
merely a detail of how the quantum tunneling is put into 
the wavefunction, but the overall picture of the resulting 
state is the same. Whether such state is energetically 
favorable compared with other competing states requires 
a more detailed study. 

Finally, we expect that the spin dynamics is highly 
degenerate in such charge frustrated systems, and its ul- 
timate fate is resolved only at much lower energy scales. 



E v /V - (E v /V) min , Jastrow-Gutzwiller ? 
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FIG. 10: This is similar to Fig. © but for doping x = 0.35 
(fermion density p = 0.65). The repulsion energy is referenced 
to the minimal repulsion energy. At this density, W c ~ 2.8. 

B. Doping x=0.35 

We now summarize similar tV study at x — 0.35. This 
is of interest for the hydrated compound Nao.3sCo02 • 
I.3H2O that was found to exhibit superconductivity. 

Figure ^| shows the expectation values of the kinetic 
and nearest neighbor repulsion energies in the ^ jg eval- 
uated using VMC. The repulsion energy is referenced to 
the minimal repulsion energy at this density [E nn rnin = 
V(3p— 1) per site; cf. Fig.0]. The renormalized meanficld 
approximation is also shown and is fairly accurate. 

The result of the wavefunction optimization for the tV 
Hamiltonian is shown in Fig. ^2 For V/t < 4.2, the op- 
timal W op t remains < 1.5; the wavefunction is metallic 
with weak renormalizations. For larger V/t, the opti- 
mal 'FjG jumps to the a/3 x y3 ordered side; however, 
the optimal W opt remains fairly close to W c w 2.8, and 
the system retains significant part of the original kinetic 
energy. 

Turning to the regime of large V, we consider also the 
more familiar CDW trial state obtained from the mean- 
field Hamiltonian Eq. (JHJ . The optimized energetics with 
such conventional CDW wavefunction is also shown in 
Fig. m The optimal Acdw remains close to zero for 

V < 3, but becomes significant and negative for larger 
V . We now observe that in the regime with the pu- 
tative a/3 x a/3 charge order both the JG wavefunction 
and the meanfield CDW wavefunction give very close op- 
timal energies. This is also true for the individual t and 

V parts, suggesting that the Jastrow-Gutzwiller and con- 
ventional CDW wavefunctions give in fact essentially the 
same physical state. 

This can be understood by examining the meanfield 
CDW state. For Acdw < the B and C sublattices are 
preferentially occupied, while the A sublattice is pref- 
erentially empty. In this case, the lower two meanfield 
bands retain much of the original bandwidth even in 
the limit of large Acdw- This is because the B and C 
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FIG. 11: Optimization of the tV Hamiltonian over Jastrow- 
Gutzwiller and conventional CDW trial wavefunctions at 
x — 0.35 (cf. Fig. [7J. We also show the result for Jastrow- 
Gutzwiller honeycomb hopping ansatz. 



sites form a connected honeycomb lattice, and for large 
Acdw the two bands correspond essentially to hopping 
on this lattice. The physical state is now obtained by 
the Gutzwiller projection of this free fermion honeycomb 
lattice hopping state. But this is also roughly the picture 
of the Jastrow- Gutzwiller wavefunction in the -\/3 x \/3 
regime for this density. 

The above suggests that we also try projecting honey- 
comb hopping ansatz, since it better utilizes the B — C 
hops. However, for the range of V/t studied here, the 
uniform triangular hopping ansatz performs better, pri- 
marily since it manages to retain some of the A — B and 
A — C hopping energy. This completes our exploration 
of the y/3 x V3 order. 

Finally, we note that at this high fermion density p = 
0.65, unlike the case with p = 0.30, the spin degeneracy 
does not occur, and the unpolarized wavefunctions are 
always better. 

To summarize, the local energetics of the tV model 
at x = 0.35 is well captured by either the renormalized 
Fermi liquid, or the V3 x \/2> charge ordered state, de- 
pending on the value of V/t. The V3 x -s/3 state also 
has mobile fermions occupying primarily the honeycomb 
sublattice (of the original triangular lattice); however, 
since the fermion density is close to complete covering of 
the honeycomb lattice, the fermion hopping is strongly 
suppressed. 

Again, the ultimate fate of the Fermi liquid (or the liq- 
uid part in the y/3 x -\/3 regime) is resolved only at lower 
energies. In the following section we study the supercon- 
ducting instability due to the J term, and whether the 
RVB superconductivity can be significantly enhanced by 
the discussed strong kinetic energy suppression. 



VI. RVB SUPERCONDUCTIVITY: 
RESURRECTION NEAR as = 1/3? 

We now turn to the issue of RVB superconductivity 
due to the antiferromagnetic spin interaction at dopings 
< x < 0.4. In the context of the triangular lattice t J 
model, this was considered by several authorsii£ii£ii2i2£i2£ 
These studies predict d + id superconductivity. As ex- 
pected for such RVB scenario, the superconductivity is 
strongest (relative to the metallic state) near half filling 
x = 0, where the charge mobility is very low. Away 
from half-filling at moderate dopings, the need to satisfy 
the kinetic energy of the carriers leads to strong suppres- 
sion of the superconductivity. As pointed out in Ref. 
and discussed further below, the experimentally observed 
superconductivity at doping x — 0.35 represents a sig- 
nificant problem to this scenario: If one uses the LDA 
bandwidth to estimate |£bare| ~ 50 — 100 meV, and takes 
the hopping integral sign as in this work, and makes a 
reasonable guess J ~ 10 to 20 meV, the resulting RVB 
superconductivity is vanishingly weak for this doping and 
would not be observed. 

As discussed above — cf. Fig. [3] — charge frustration can 
lead to strong suppression of the effective hopping ampli- 
tude i e ffi even for larger doping. Here we study whether 
this suppression can be strong enough to resurrect the 
superconductivity at x = 0.35. Figure El also suggests 
that the region near x = 1/3 is special in that it allows 
the strongest such renormalizations, with or without the 
charge ordering. As discussed earlier, this is because the 
charge system is most sensitive to the further neighbor in- 
teractions near this commensurate doping. On the other 
hand, when the superconductivity is weak, the transi- 
tion temperature is exponentially sensitive to the effec- 
tive hopping amplitude (see below). Thus, we may spec- 
ulate about the possibility of a small superconducting 
dome around this special doping due to charge correla- 
tion (possibly, charge ordering). 

The physics treatment presented below is very 
schematic. We will essentially think only in terms of 
the renormalized couplings and ignore the fact that the 
underlying "liquid" state may be charge-ordered. This 
is done to get a rough feeling as to whether the sug- 
gested scenario can work at all. If the superconductivity 
near x = 1/3 can reappear only in the charge-ordered re- 
gion, the obtained insight is still useful and may suggest 
a more careful treatment. For example, the charge order 
may be suppressed by longer range repulsion. Another 
possibility is that a more accurate treatment may lead to 
short range charge order for intermediate coupling which 
shares the same kinetic energy suppression. 

We first summarize the standard RVB meanfield for 
the pure t J model. We formulate this meanfield as an 
approximate variational procedurei2ii22i2& This is partic- 
ularly convenient for the present work, which also takes 
the variational wavefunction perspective. 

To study the possibility of singlet superconductivity, 
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we consider "trial" Hamiltonian 

Atrial = [-Xijcl.Cja - (Ayc]^ +h.C.) 

ij 

- X ^ ~ ^ ~ ^ ' 

with = x*ji = Aij. For each such trial Hamilto- 
nian, we obtain the corresponding ground state |^o)- In 
the meanfield, we ignore the no-double-occupancy con- 

I 



straint and only require the average density to be cor- 
rectly {(y i(J Ci a ) = 1 — x, which is achieved by tuning the 
chemical potential fi. Going beyond the meanfield, the 
physical wavefunction is obtained by Gutzwiller projec- 
tion. 

As discussed at length earlier, we can approximate the 
expectation value of the t J Hamiltonian in the physical 
wavefunction by proper renormalizations of the meanfield 
values: 



QE G |g fJ |*< 



>9t- 



(\ Eolff«lg ) 
(*o|*o) 



-gj- 



(tt |ffj|*o) 



(*ol*o) 



-9t^2Uj{c\ a c ja ) +c.c. 

(ij) 



-gj 



(ij) 8 



\{c\ a c ja )\ 2 + \{e aa ,c]J 



(10) 



The hopping renormalization factor is given by Eq. Ipjjl. 
while for the Heisenberg exchange we haveSLS 2 . 



These estimates follow essentially from the no-doublc- 
occupancy configuration constraints, and do not depend 
on the details of the preprojected state as long as it is 
spatially uniform. Also, they give numerical results that 
are fairly close to the actual evaluations with the pro- 
jected wavefunctions, as discussed earlier. The above 
is pr ecisely th e renormalized meanfield formulation of 
Refs. I2ll22l26l The slave boson meanfield of Ref . |l9| uses 
instead gt — x and gj — 1, so their numerical values are 
somewhat different. 

In this formulation, only the ratio A/\ is meaning- 
ful. A convenient procedure to minimize Eq. ||Tg is to 
minimize instead the so called meanfield Hamiltonian 

^ mf = X q — r [IXtt ~ 9tt tJ \ 2 + I Ay | 2 ] + Atrial . 

(ij) 69sjij 

By standard arguments, the global minimum of the 
meanfield Hamiltonian is also the minimum of the trial 
expectation value Eq. I|10l) . In this formulation, the op- 
timal x an d A each obtain physical scale as set by t 
and J. Thus, we can get a rough idea about the quasi- 
particle spectrum above the ground state by considering 
the meanfield excitation spectrum, which now has proper 
physical scale. In particular, the optimal A gives a phys- 
ical measure of the strength of superconductivity. 
The self-consistency conditions read 

X% = 94 + ^f(cU,), (12) 
A*, = • (13) 

From now on, we specialize to the d+id superconductor 



ansatz: 

A ei = A; A e2 = Ae' 21 / 3 ; A e3 = Ae l4 ^ 3 . (14) 

Here e\ — x, e 2 — \x + ^y, and e 3 = e 2 — e\ 
are the unit triangular lattice vectors. There is strong 
evidence that this state wins the t J model energetics 
for the considered dopings, at least on the meanfield 

lrT l 17.18.19.26.27 

We give the results for t = 3J and t = 5J. This is 
somewhat different from the cited t = 5 — 10 J values^ 
At the moment, there is significant uncertainty in the 
precise microscopic model, while the superconductivity 
energy scale is exponentially sensitive to the microscopic 
values and to numerical constants in the theory. The 
t = 3 J results make the demonstration of principle more 
dramatic. A similar, but weaker, effect is seen for i = 5 J 

The optimal A in units of J is shown as a function of 
doping in Fig. 1121 For weak J, the optimal A is given by 
a BCS-like formula (see Appendix lAl for details) 

A ~ t cS e- ct ^ /J ^ (15) 

with some numerical constant c = c(x) and t c g = g t t, 
Jeff — gjJ- The effective mobility of charges increases 
with doping t c g ~ xt, and this leads to the observed very 
quick drop of A. 

Figure^Jalso shows two other physical measures of the 
strength of superconductivity. One is the meanfield T c 
defined here as the transition temperature for the finite 
temperature optimization of the meanfield Hamiltonian 

The other measure is obtained by considering the con- 
densation energy of the superconducting state. This is 
defined as the energy gain in the optimal superconduct- 
ing state relative to the Fermi liquid state (A = 0). For 
small A, the condensation energy is expected to scale as 
-Econd ~ A 2 /t fj; to compare with A in Fig. [21 we plot 
instead (-Bcond^/) 1 ^ 2 • 
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FIG. 12: Renormalized meanfield results for the d + id super- 
conducting state for t — 3J. We show the meanfield T c , the 
optimal A, and the square root of the condensation energy 
Scond- The energy scale is the bare J; note the logarithmic 
scale for the energy. 
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FIG. 13: Jastrow-Gutzwiller renormalization factor for the 
Heisenberg exchange as a function of doping for a number of 
fixed W (cf Figs. I3I4H . In the \/3 x \/3 phase (data points 
below the thick dark line), the exhibited gj is averaged over 
all bonds. 



From Fig. 1121 these measures all trail each other. The 
figure has been somewhat arbitrarily cutoff at 10~ 3 J. For 
J ~ 20meV, any T c or A below this scale would not be 
observed in the experiments. Note the precipitous drop 
in the strength of superconductivity for x > 0.20. There 
is simply no hope for it surviving to the experimentally 
observed x = 0.35 in this setting. 

We remark here that a direct VMC study must see 
the condensation energy to establish the ground state 
A. Since E con d is extremely small, such studies become 
impractical. This is where the renormalized meanfield 
procedure becomes very useful. 

Let us now return to the tJV model with strong near- 
est neighbor repulsion. We think roughly as follows. The 
dominant t and V parts can be satisfied as above by the 
appropriate Jastrow weighting of charge configurations in 
our trial wavefunctions. As discussed earlier, the effect of 
the Jastrow factor can be conveniently described by the 
corresponding renormalizations of the hopping amplitude 
9t (Eq. and Fig. and the Heisenberg exchange gj. 
The latter is approximated by 
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and is plotted in Fig. [23 (cf. discussion following Eq. 0. 
As long as the charge distribution remains uniform, these 
renormalizations capture the main effect of the nearest 
neighbor correlations built in by the Jastrow factor. 

Thus, for each doping level x and the Jastrow sup- 
pression strength W, we can estimate the correspond- 
ing t e s, Jeff, and then the optimal A. The latter is our 
main measure of the superconductivity strength and is 
shown in Fig. The W = line is the same as in 
Fig. 1121 while the W = 8 corresponds essentially to the 
minimum-nearest-neighbor projection. Again, the dark 
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FIG. 14: Renormalized meanfield for the Jastrow-weighted 
d + id superconducting state for t — 3J (cf. Fig. 11211 . We show 
the self-consistent A as the measure of the superconductiv- 
ity strength (T c plots look very similar). The dark thick line 
corresponds to the phase boundary of the Jastrow weight (cf. 
Figs. I2I3I13H — the maximum enhancement of the supercon- 
ductivity while remaining in the uniform phase. 



thick line corresponds to the phase boundary of the Jas- 
trow weight. For x > 0.27 all points above this line have 
the v3 x y3 charge order. These are obtained by using 
the corresponding formal renormalization factors and the 
above prescription, even though this violates the initial 
motivation coming from a uniform renormalized liquid 
picture. As emphasized earlier, the precise energetics 
in this regime likely requires a more careful treatment. 
However, we expect that even such simplistic analysis in 
the x V3 regime gives a reasonable first guidance on 
how the role of the kinetic energy can be suppressed by 
possible charge ordering in the system. 

Our tentative conclusion from Fig.^Jis that for t = 3 J 
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FIG. 15: This is the same as Fief. ITU but for t = 5 J. 



view the fermions as restricted primarily to the honey- 
comb lattice. On the other hand, the d + id state wins 
the energetics in the original uniform triangular lattice 
meanfield and needs to be re-examined in the present 
context. The above renormalized meanfield procedure 
roughly corresponds to restricting the d + id ansatz onto 
the honeycomb lattice. Of course, one should also con- 
sider other possible RVB superconductor states on the 
honeycomb lattice and decide which one is optimal ener- 
getically. More generally, one may want to consider trian- 
gular lattice superconducting ansatz with broken transla- 
tional symmetry patterned after the -s/3 x state. We 
are not pursuing such studies here, since it is important 
to first establish whether the charge ordering occurs at all 
in the material. If this indeed happens, the above rough 
considerations can give us some initial idea about the 
scale of the superconducting instabilities in such state. 



the considered nearest neighbor Jastrow renormalizations 
that leave the underlying liquid wavefunction in the pure 
Fermi liquid state are borderline sufficient to explain the 
superconductivity near x = 1/3. One should of course 
judge this critically because of the exponential sensitivity 
to the actual value of the ratio t e s/J e g, Eq. I)15[l. The 
trend for increasing t/ J can be seen by comparing t = 3 J 
Fig. [Hand t = 5 J Fig. [H 

We speculate that in Naa;Co02 • I.3H2O the actual sit- 
uation is qualitatively close to the curve with W = 3 in 
Fig El which near x = 1/3 roughly corresponds to the 
critical strongly correlated liquid of the nearest neighbor 
Jastrow weight (cf. Figs. I2I3J) . This curve may lie inside 
the liquid phase for longer ranged Jastrow weight. An- 
other possibility is that only short range charge order is 
developed for intermediate coupling. One thing should 
be clear from Figs. 1141151 there can indeed be signif- 
icant enhancement — one to two orders of magnitude — 
in the superconductivity scale due to the kinetic energy 
suppression in the charge-correlated liquid. Because the 
charge system is most responsive near the commensurate 
x = 1/3, this enhancement may be strongest near this 
doping, which may explain the experimentally suggested^ 
superconducting "dome" around x = 1/3. However, we 
note that even for W = 3, A in Fig. [21 shows only a 
shallow maximum near x = 1/3. A possible explanation 
of the experiment is that x significantly less than 1/3 is 
not achievable due to chemical reasons and superconduc- 
tivity is simply cut-off. Again, the important message we 
draw from Fig. is the possibility of pushing T c up to 
an observable level near x — 1/3. 

Finally, if we continue the theory into the Jastrow 
V3 x \/3 charge order regime, the T c enhancement may 
be even stronger. This is not surprising, since the charge 
mobility is suppressed even further in this case. Thus, 
our earlier analysis tells us that for x > 1/3 we are es- 
sentially doping a nearly half-filled honeycomb lattice. 
This picture also suggests some possibilities of treating 
the y/3 x y/3 regime more carefully, similar to our dis- 
cussion in Sec. [V] For example, for x > 1/3, we can 



VII. CONCLUSIONS. CONNECTION WITH 
EXPERIMENTS. 

We conclude by stating some consequences of the dis- 
cussed effects of charge frustration. 

1) It will clearly be interesting to look for signs of 
charge order near x = 1/3 and 2/3 using X-ray or neu- 
tron scattering. The conductivity is metallic and in the 
case of x — 1/3 reaches 50e 2 /h at low temperatures^ 
This suggests that long-range charge ordering is unlikely, 
but there may be a tendency for short-range ordering. 

2) There can be strong suppression of the effective hop- 
ping amplitude due to nearest neighbor repulsion while 
remaining in the Fermi liquid state. The meanfield hop- 
ping amplitude Xij = X is given by Eq. (|12fl and has 
contributions proportional to gtt and gjJ. Note that in 
addition to the suppression of gt (Fig. [3, gj is also sup- 
pressed (Fig. IT3|) . especially for x > 0.5. 

This suppression leads to low fermion degeneracy tem- 
perature. The properties of such Fermi liquid system 
with tp <T are rather unusual from the perspective of 
the familiar metals with ej? » T (the Fermi energy is 
measured from the bottom of the band, and is roughly 
£F ~ t Q g). 

a) In particular, the thcrmopower is large and satu- 
rates to the value 



jx_ _ ks_ 2 - p 



qT 



9 



(17) 



at large temperature. Note the "classical" scale fes/|e| = 
86.2 nV/K, which is in fact observed in Na 2 .Co02 
The full temperature dependence for x = 0.70 is shown in 
Fig. lib! Here and below, we use simple-minded transport 
theory summarized in Appendix [5] Form Fig. 1161 the 
thermopower reaches one half of the maximal value for 
T as t e ff. 

b) The Hall coefficient for the triangular lattice band 
structure has an unusual non-saturating increase with 
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FIG. 16: Simple-minded transport theory for nondegener- 
ate spin-1/2 Fermi gas on the triangular lattice. The ther- 
mopower is plotted in "classical" units of fcs/|e|, while the 
units for the Hall coefficient contain the three-dimensional 
volume per Co atom fl. 



the temperature for T > tp as observed experimentally 
in Ref. llOt The limiting high temperature behavior is 



Rh 



qcp(2-p) t cfi 



(18) 



where fi is the three-dimensional volume per Co. The 
full temperature dependence is shown in Fig. 1 161 the high 
temperature trend sets in for Twlto 2 t e ff. Possibility 
of this unusual behavior was predicted in Ref. from the 
high temperature expansion for the t J model (the doping 
dependence of the proportionality coefficient is somewhat 
different here). We remark that this unusual behavior is 
the consequence of the triangular lattice band structure 
only, and its origin can be traced to the presence of three- 
site hopping loops as detailed in Appendix 151 Correlation 
effects per se are needed only to reduce t g below the 
experimental temperatures. 

c) Pauli susceptibility per Co site for T > ep becomes 



model (related to the lattice gas with nearest neighbor 
repulsion as mentioned in Sec. Illlfl . Note that the band- 
width observed by heat capacity^, and by ARPESii is 
proportional to x an d has contributions from both g t t 
and gjJ (see Eq. I12jl , On the other hand, electromag- 
netic response couples only to t so that the Drude weight 
observable from infrared reflectivity and the superfluid 
density (observable via the London penetration depth in 
the case of superconductors) are directly proportional to 
gt- These will provide a more sensitive test of the pre- 
dicted dip in g t near x = 1/3 and x = 2/3 as shown in 
Fig. |21 For example, it will be interesting to compare the 
Drude weight for x = 1/3 and x = 0.5 samples. 

4) The spin physics near x = 2/3 is expected to be 
highly degenerate and complicated, and will manifest it- 
self below the energy scale t c g. In particular, the above 
transport pictures will likely be modified below this scale. 

5) Near x — 1/3, whether the system prefers uniform 
or charge-ordered state, the correlated liquid can have 
further RVB superconducting instabilities. The suppres- 
sion of the charge mobility serves to enhance and may 
even resurrect the superconductivity under a small su- 
perconducting dome around x = 1/3. Experiments^ ob- 
serve the disappearance of the superconductivity below 
x = 0.26. However, the strongest RVB superconductivity 
is expected at much lower doping, and the search should 
be pursued more vigorously towards x = 0, if that is 
chemically possible. 
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(19) 



Note that this has a Curie-like behavior, but is some- 
what smaller — by a factor of 1 — p/2 — from the case of 
completely free spins. 

3) The kinetic energy renormalizations are strongest 
(for fixed repulsion strength) near the commensurate 
x = 1/3, 2/3, and weakest near x = 1/2. This is be- 
cause the system finds it easiest to order, even if only 
locally, near the commensurate filling, while away from 
commensuration much of the nearest neighbor repulsion 
energy cannot be avoided in any case. 

Charge frustration may also be relevant for the ex- 
perimental "charging" curve of Ref. 0. The observed 
plateaus at x = 1/3, 2/3 remind one of the magneti- 
zation plateaus in the frustrated triangular lattice Ising 



APPENDIX A: DETAILS OF EQ. lfl5|) 

Eq. (|15fl can be understood by examining the self- 
consistency conditions Eq. (|13|) . Specializing for the d+id 
ansatz, we have 



1 = 



3J C 



ff 1 \ - 



fd+id{k) 



fd+id{k) = 2 cos k ■ ei ( cos k ■ e\ — — cos k ■ e 2 — — cos k ■ e 3 



Here Ni at t is the number of lattice sites; = — p with 
£fc = — 2%(cos k ■ ei + cos k ■ e 2 + cos k ■ e^). 

For weak superconductivity A <^ J c ff < £ ff , following 
a BCS-like analysis, we obtain the following approximate 
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Triangular lattice: ingradients for A(0), T c 




FIG. 17: Density of states per site vo(E) and the integrated 
DOS N(E) for the triangular lattice band structure (unit hop- 
ping amplitude, no spin). We also show the data for the aver- 
age d + id wave factor fd+id at the corresponding energy cut; 
we plot 0.5/d+id to fit into the same vertical scale. To ob- 
tain the values corresponding to the particular doping x, we 
first locate E such that N(E) = (1 — x)/2. In this manner, 
we obtain c(x) — 4/(3^0 /d+id) plotted for the relevant range 
< x < 0.4 in the inset. 



a) The thermopower is given by 

Q = — ■ (Bl) 

&XX 

The kinetic coefficients are given by the integrals over the 
Brillouin zone 

[ d 3 k ( df\ ,,,e(fc)-,u 

2 r d 3 k ( df\ „, „, 

^ = « T _/ ^3 J «x(*)«x(*) • 

Here, /(e) = l/(e e_A1 + 1) is the Fermi distribution. 
We cite the more familiar three-dimensional expressions. 
When applying to Na^CoC^ , we specialize to the layered 
triangular lattice by assuming no dispersion in the z di- 
rection. The full temperature dependence for x — 0.70 
is shown in Fig. 1161 and the limiting high-temperature 
behavior is given in Eq. 117|l . 

b) In weak magnetic fields uj c t -C 1, the Hall coefficient 
is given by 



Rh = 



OH 



' xx w yy 



(B2) 



formula 



A = A x exp 



?,V (x)fd+id{x) JeS 



(Al) 



A is an order one numerical constant, frj(x) is the trian- 
gular lattice hopping density of states per site (not in- 
cluding spin) at the Fermi energy corresponding to dop- 
ing x, and fd+id is the d + id wave factor averaged over 
the Fermi surface. The scale x ln front of the exponen- 
tial corresponds to the energy cutoff being roughly the 
Fermi energy, since the pairing is over the full Fermi vol- 
ume. For small J e s we see from Eq. (|12|) that x can 
be replaced by t cff in (|A1I) . yielding Eq. (f[3)l . Similar 
expression is obtained for the meanfield T c . 

The necessary data is shown in Fig. El The coefficient 
c(x) depends rather weakly on x, and the main effect 
on the A and T c is from the doping dependence of t e s- 
The above approximate formula agrees fairly well with 
the actual meanfield calculations performed in the main 
text. 



APPENDIX B: TRANSPORT FOR T < e F 

In this appendix, we summarize the simple Fermi liquid 
transport theory that was used to obtain Fig. ^| and 
Eas. ll7lf81 The main formulae can be found in standard 
texts. 28 - 29 . 



o X x , Oyy are the static zero field conductivities given ear- 
lier, while 



q 3 T 2 



In the last equation, M Q 1 (fc) = 
mass tensor. 



dk a dkfj 



(B3) 
is the inverse 



The high temperature behavior for the layered trian- 
gular lattice is given by Eq. I|18f) . The origin of this non- 
saturating increase with temperature lies in the presence 
of triangular hopping loops. Indeed, consider the above 
semi-classical expression for an at high temperature, and 
translate it from the momentum space back to the real 
space assuming a general hopping problem t&R/ on a Bra- 
vais lattice. The result reads 



Off 



Here, 



q 3 T 2 



df 

de 



de 



E 



^01^12^20 RlxR2y{Rl><R2 



p(2 — p)f (4T) , and also enters a xx , a y 



(B4) 
; ft is 



the volume of the unit cell. The lattice hopping problem 
is input through the real space sum over possible hops out 
of the origin: Ri = i?oi , Ri = R02 ■ For each triangle 
specified by an unordered triple of vertices 0, R\, R2, the 
clockwise — > 1 — > 2 — * and anticlockwise 0^2^ 
1 — > contributions add to (i?i x R^H, i.e., a quantity 
of definite sign. The effect is of course strongest for the 
triangular lattice. 
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